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Abstract 



It is known that Newton's law of gravity holds asymptotically on a flat "brane" embedded in an anti-de Sitter "bulk" ; 
this was shown not only when gravity in the bulk is described by Einstein's theory but also in Einstein "Lanczos 
Lovelock Gauss-Bonnet" 's theory. We give here the expressions for the corrections to Newton's potential in both 
theories, in analytic form and valid for all distances. We find that in Einstein's theory the transition from the I/r 
behaviour at small r to the 1/r 2 one at large r is quite slow. In the Einstein Gauss-Bonnet case on the other hand, 
we find that the correction to Newton's potential can be small for all r. Hence, Einstein Gauss-Bonnet equations in 
the bulk (rather than simply Einstein's) induce on the brane a better approximation to Newton's law. 

I. INTRODUCTION 

There has been of late an increasing interest for gravity theories formulated within spacetimes with one "large" 
extra dimension, that is, the idea that our universe may be a four dimensional singular hypersurface, or "brane" , 
embedded in a five dimensional spacetime, or "bulk". The Randall-Sundrum scenario p], where our universe is a 
four dimensional quasi-Minkowskian edge of a double-sided perturbed anti-de Sitter spacetime which satisfies the five 
dimensional Einstein equations, was the first explicit model where the linearized four dimensional Einstein equations 
were found to hold on the brane, at least far from the sources, Newton's potential being corrected by a small 1/r 2 
term. This claim was thoroughly analyzed and the corrections to Newton's law exactly calculated, either in the small 
r limit where it was shown that they diverged in 1 jr or in the large distance limit . The analysis was then extended 
to the case when gravity in the bulk is described by the Einstein Gauss-Bonnet equations, and it was found that, 
in that theory as well, Newton's law also held on the brane at large distance, 0, Q with various results on the exact 
1/r 2 corrections, due to more or less careful treatments of "brane bending" and conflicting views on the boundary 
conditions to be imposed on the brane, see e.g. for a review. 

In this paper we describe gravity in the five dimensional bulk by means of the Einstein Gauss-Bonnet equations 
and use the junction conditions obtained in [(| (which are briefly reviewed in section 2). To obtain the geometry of an 
almost flat brane imbedded in a perturbed anti-de Sitter bulk and, hence, its gravitational field, we use the approach 
advocated in 0. 

More precisely, we first construct in section 3 the most general linear perturbations of the bulk satisfying the 
Einstein Gauss-Bonnet equations ; it turns out, as already known, 0,0 that they are the same as in Einstein's theory 
(the bulk background being locally anti-de Sitter). Then, in section 4, we use the junction conditions obtained in 
[6j and write the equations which govern gravity on the brane ; we deduce from them a set of equations that can 
be compared to the usual four dimensional linearized Einstein equations ; we see, as already known, [3J, |^| that the 
two sets (after proper identification of Newton's constant)are identical if only "zero modes" are present in the bulk 
(however only very special matter sources on the brane are compatible with zero modes only). 

We then turn, in section 5, to the case when matter on the brane is a static, spherically symmetric, point source 
andgive the linearized gravitational potentials in terms of well defined integrals which generalize the results obtained 
in pj (equations H5.6l5.8f) in text). We finally present our main results in sections 6 and 7 : by approximating the 
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integrands by simpler functions we are able to perform the integrations and thus obtain the gravitational potentials 
in analytic form valid for all distances. 

This allows us to see that, when Einstein's theory governs the bulk, the "correction" to Newton's potential on the 
brane goes from its large 1 jr behaviour at small r to its 1 jr 2 one at large r quite slowly ; in the Einstein Gauss-Bonnet 
case on the other hand, we find that the correction to Newton's potential can be small for all r. We hence conclude 
that the Einstein "Lanczos Lovelock Gauss-Bonnet" equations in the bulk (rather than simply Einstein's) induce on 
the brane a better approximation to Newton's law and the usual linearized /o«r-dimensional Einstein equations. 

Throughout this paper, we denote double partial derivatives d^dy by d^ v . 



II. THE EINSTEIN "LANCZOS LOVELOCK GAUSS-BONNET" EQUATIONS AND THE 

ASSOCIATED JUNCTION CONDITIONS 

To construct a "Z 2 -symmetric braneworld" one can proceed as follows : consider a five dimensional Lorentzian 
manifold V+ with a timelike edge, and superpose V+ and its copy V- onto each other along the edge ; one thus obtains 
a braneworld, composed of a ^-symmetric spacetime or "bulk" V5, and a singular surface, or "brane" £4 whose 
extrinsic curvature is discontinuous : the extrinsic curvature of £4 embedded in V— is the opposite of the extrinsic 
curvature of £4 embedded in V+. Some components of the braneworld Riemann tensor therefore exhibit a delta- type 
singularity in V5 U £4. In a Gaussian normal coordinate system such that the equation defining the position of the 
brane is y — the braneworld line element can be written as 

ds% = dy 2 +g^(x p , \y\)dx p dx v (2.1) 

where y > spans V + and y < spans V— . As for the brane extrinsic curvatures in V + and V_ , they are given by 

= K+ u = -K~ v = -i%| y=0+ . (2.2) 



2 dy 



We shall associate to this braneworld the following action 



S= / d 5 xy^g~L [2] +2K d 4 xy^§C m -2 d 4 xy^jQ. (2.3) 

JVs JS 4 JY.A 

g is the determinant of the bulk metric coefficients qab , 9 that of the induced brane metric coefficients g^ v . In the 
first term : 

L [2 ] = -2A + s + aL (2) , 

where L (2) = s 2 - Ar LM r LM + R LMNP Rlmnp , (2.4) 

is the Einstein Lanczos Lagrangian (also called Einstein Gauss-Bonnet and generalized by Lovelock ||, see e.g. 
for a review). Here R A bcd = dc^go — ■■■ are the components of the Riemann tensor, T^ D being the Christoffcl 
symbols and all indices A being moved with g&B and its inverse g AB ; tbd = R A bad are the Ricci tensor components, 
s = g BD rBD is the scalar curvature ; A is the bulk "cosmological constant" and a a (length) 2 parameter (that we 
shall take to be positive, see e.g. 9]). In the second term, C m {g^ v ) is the brane "tension plus matter" Lagrangian 
and k a coupling constant which will be ultimately related to Newton's. The third, boundary, term, which generalizes 
Gibbons-Hawking's was ^ rs ^ obtained by Myers and its explicit expression is (see Q) 

Q = 2K + 4a(J - 2a£ K£) (2.5) 

with J the trace of 

J£ = -Ik$KPKZ + \KK$K$ + l -K»{K.K - K 2 ) ; (2.6) 

<j£ = fj} — \5vS is the intrinsic Einstein tensor of the brane £4 and Kjf its extrinsic curvature (|2.2|l . all indices \i 
being moved with g^ v and its inverse g^ v . 

Thanks to this boundary term the variation of S with respect to the metric coefficients is given in terms of their 
variations only [llj, as 0] : 

SS= f d b x^—go [2 } B Sg AB + f d 4 x^§ (2fl M „ - kT^) 5g^ . (2.7) 

JV 5 J £4 



3 



The "braneworld equations of motion" are therefore simply taken to be 5S = 0, with the metric fixed at the boundaries 
at infinity only. They are, first, the Einstein Gauss-Bonnet "bulk" equations (which are second order and quasi-linear 
in the metric coefficients, see e.g. p| for a review) : 

^]B = ^B + ^ + aaf 2)B =0, (2.8) 
where erg = r B — \5 B s is the bulk Einstein tensor and 



a 



A )B = 2 [R ALMN R BLMN - 2r LM R A LBM 2r AL r BL + sr A ] U A L (2) (2.9) 



is the Lanczos tensor. As for the brane equations, which, as shown in 6|, generalize the Israel junction conditions [T^ 
to Einstein Gauss-Bonnet gravity, they are : 

Bj} = KZ- K5» + 4a (~J<f - IjSfJ - P^kA = fa (2.10) 



where 



Pp,pvcr — Rf-ipva H~ (j'fKjdpi' ^pcrSfif H - ¥ ' pvQpia ¥ pyQ per) 

~2 S (9»<t9 P u ~ 9nv9 P o) (2-11) 



and where T^ v is defined by S(y/^g£ m ) = ~\\/—gT^ v bg^ v and is interpreted as the stress-energy tensor of "tension 
plus matter" on the br 
()2.1(l . one can see that 



2 ' 

plus matter" on the brane. By making a (4 + 1) decomposition of <J^ B in the Gaussian normal coordinate system 



tfq* = ( 2 - 12 ) 
with V;, the covariant derivative associated with the brane metric g^, so that the conservation of T^, 

V !/ T^ = 0, (2.13) 



is included in the bulk equations 

If now the bulk is locally an anti-de Sitter (AdS5) spacetime, then, because of maximal symmetry : 

Rabcd = --^{gAcgBD - gADgBc) , (2-14) 

with 



1 1 / / 4aA \ 

^=4^( 1± V 1 + — J (2 - 15) 

in order to satisfy the bulk equations Ij2.8|l . (One usually chooses the lower sign so that lim a ^o C 2 = -f , that is the 
Einsteinian value ; however, there is no reason at that level to reject the upper sign which allows, in particular, for 
an AdSs bulk even when A = ; as for the case C 2 = 4a it has been analyzed in e.g. H3-) moreover, one wants 
the brane to be flat, a convenient coordinate system to describe the anti-de Sitter bulk is 

ds% = dy 2 + e- 2 ^/ c Vll „ dx^dx" (2.16) 

in which the brane is located at y = and has an extrinsic curvature given Ktf — see (I2.2|) . The brane equations 
(|2~TU|) then tell us that 

6 / 4a \ 

^--c^-w) (2 - 17) 

which is interpreted as the "tension" required for keeping a brane flat in an anti de-Sitter bulk governed by the 
Einstein Gauss-Bonnet equations. 
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III. THE BULK PERTURBATIONS 



In (quasi) conformally Minkowskian coordinates x A — {x p = (t, x l ),x i — w}, the line element of our perturbed 
.^-symmetric anti-de Sitter bulk can be taken to read 

ds 2 \ 5 = f^j (tjab +jAB)dx A dx B with j ww = 7^ = (3.1) 

where C{> 0) is given by l|2.15[l and where 7^„ are ten fonctions of the five coordinates (x M , w). Since the (unperturbed) 
brane is located at w = 0, the coefficients 7^ will be at first order functions of \w\ only. 

If we now impose this metric to satisfy the bulk Einstein Gauss-Bonnet equations (|2.8|) . it is an easy exercise to 
see, first, that the conditions (all indices being from now on raised with r] pv ) 

7 P P = , d„7£ = (3.2) 

solve the constraints &fy w = 0, c^]^ = at linear order and reduce the ten metric coefficients 7^ to five, which 
represent the five degrees of freedom of AdSs gravitational waves in Einstein Gauss-Bonnet theory (this being so 
because the Lanczos-Lovelock Lagrangian yields second order field equations, just like Einstein-Hilbert's). 
As for their evolution equation, it then reads (with Qi = ii pu d p d v ) 

1 - ^ j ( D47m^ + dwwl^u - ^-dwl^v J = . (3.3) 



w 



The presence of the overall coefficient (1 — 4a/ C 2 ) = ^fy/T+ aA/3 is a consequence of the quasi- linearity of the 
Lanczos-Lovelock equations <|2.8[) . see e.g @- The case when the parameters a and A are such that it vanishes 
(studied in e.g. ,f{|) will not be considered here. Therefore, the evolution equation for the AdSs gravitational waves 
in Einstein and Einstein Gauss-Bonnet theory are the same 0, Q . 

As emphasized in 0, the general solution of (|3.3|) must converge (more precisely be L2) on its domain of definition. 
Now, the bulk is composed of V+ and its copy V- joined by the brane £4 and its Z2— symmetric metric (|3.I|I depends 
only on \w\. The general solution of l|3. 3fl is therefore the union of two sets. First, the set : 

ltlv {x»,w) = Tie /-^-^re^'e^^H^I) (3.4) 
J (2tt)2 (2tt)2 



with k = -\/k 2 +m 2 , k p e pfl = , r\ p(I e prj = 

where k 2 = k, L k l , where the five polarisations e pu are a priori arbitrary functions of k l and m and where Z2(mw) = 

H2 (raw) + a rn H2 2 \mw) is an a priori arbitrary linear combination of second order Hankel functions of first and 
second kinds |14|. The coefficient a m is determined by the model at hand and one usually eliminates the mode 
coming from +00, that is sets a m = . The "zero modes" do not depend on w and are : 

J (2ir)2 

with ke 0fJ _ + k l e lp = , ff° 'e pa = . (3.5) 
The second set is composed of the extra L2 modes which converge exponentially as \w\ — > +00 : 

lpv {x»,w) = nej-^j f -^re^'V^^aVPlH) ( 3 - 6 ) 

J (27r)2 J- k 2 (27r)2 

with fc = -\/k 2 + A , k p e pp = , r] pa e p<J = Q. 

The static modes do not depend on time and are : 

1 ${x\w)=ne [ ^e ik ^e p ^)w 2 H^(ik\w\). (3.7) 
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IV. THE LINEARIZED EQUATIONS FOR GRAVITY ON AN EINSTEIN GAUSS-BONNET BRANE 

The position of the brane £4 in the bulk V5 is defined in the coordinate system (|3.1|l by 

w = C + C{x») (4.1) 

where the function C,{x^) is a priori arbitrary and describes the so-called "brane-bending" effect. The induced metric 
on E4 is, at linear order 

ds 2 | 4 = {rj^ + h^ LU )dx^dx v with = 7 m „|e — 2— T)^ (4.2) 

where r Y flv (x IJ ',w) is a solution of I|3.2I3.3[) and where the index E means that the quantity is evaluated at w = L. 
As for its extrinsic curvature (defined by (|2.2|) in a Gaussian normal coordinate system) it is given, in the system 
ED, by El 

K ^ = 2 6 "~\ dwl ^ 11 + d ^- (43) 

Now, as we have seen in section 2, the "brane equations" are H2.10JI . In the case at hand Kf} is given by (|4.3(l and 
P^ pva is to be computed by means of the metric l|4.2|l . Splitting into a brane tension and matter stress-energy 
tensor S£ as 

^ = -|^(l-|)+/s5ff (4.4) 
where, from now on we shall use the notation 

4a 

a = , (4.5) 

the brane equations (|2.10|l read, at linear order 

nS£ = 2 (1 + a) (d{X - SSDiO -(!-") «|s - &CD4T£|e ■ (4.6) 

Equations l|4.2|l and (|4.6|) together with l|3.4|l and (|3.7|) completely describe gravity on the brane. They have two 
useful consequences 

d p Sg = , k3 = -6(l + a)OiC- (4.7) 

Let us compare them to the usual four dimensional linearized Einstein equations. In order to do so, we first perform 
a coordinate transformation x^ — > x^ = x*^ + e M such that the new metric coefficients 

2 2 
K-v = "V^s ~ -^CVnv + d n e » + 9 » e n witn U i e^ = -—d^C, (4.8) 



satisfy the harmonicity condition 



( h*/ - ~S£h* ) = (4-)) 



Taking the d'Alembertian of (|4.8|) . eliminating £ by means of (|4.(i|) and using I|3.2I3.3|I we get the following consequence 
of the brane equations 



1 



1- a 
1 + a 



(4.10) 



Now, recall that the usual linearized Einstein equations on a four dimensional Minkowski background read, in 
harmonic coordinates 

□4/i*„ = -16ttG - l^uS) (4.11) 
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where G is Newton's constant. One therefore sees that, if we identify, in agreement with |^| 

and if only zero modes are present in the bulk, then (|4.1U|) reduces to (|4.11ll . 

There is however an important proviso. The usual Einstein equations (|4. 1 lf> hold for any type of matter (compatible 
with the harmonicity condition, or, equivalently with the Bianchi identity d p S^ = 0. By contrast, the linearized 
equations for gravity on a brane are l|4.10|l . provided the source satisfies the junction condition l|4.6[) . If only 
zero modes are present in the bulk, then, since the last two terms in 1)4. 6J) are absent because of 1)3. 3[l . the derivatives 
dxiSftu — ^rj^S) must be symmetric in A and [i, a property which is satisfied only by very special matter, obeying 
a very contrived equation of state (see e.g. for details). There is therefore no compelling reason, at that level, to 
make the identification (4.12), that we shall consider until further notice as a mere notation. 

Now, of course, the junction conditions are better seen as boundary conditions on the allowed modes in the bulk. 
Indeed, if matter on the brane is known, then (|4.til4.7|) can be inverted to give £ and the last two terms in (|4.6|) in 
terms of S^. With that information one can in principle get, by inversion of (|3.4|) and l|3.7f) . the polarisations e^, 
and hence the allowed bulk perturbations, in terms of the brane matter source. Then the induced metric on the brane 
(14. 2|) is known in terms of the matter variables and can be compared with the usual four dimensional Einstein result. 
This programme is completed below in the particular case of a point static source. 



V. THE GRAVITATIONAL POTENTIALS OF A STATIC POINT SOURCE ON AN EINSTEIN 

GAUSS-BONNET BRANE 



Consider a static, spherically symmetric, point-like source on the brane. Its stress-energy tensor is, (with r= {x 1 }, 



Xj,X ) . 

Soo = M5(r) 
When necessary, we shall go to Fourier space : 

/ / with / 

(2tt) 3 / 2 ' r 
The junction conditions (j4.614.7l) then give : 

1 kM 1 



Soi — Sij — . 



(5.1) 



d 3 k 
{2nf/ 2 
2 1 



nk i and 



2 k 'i k j 



7T k 2 



c = - 



247r (1 + a) r 



(l-a)d wl ifk+a£A^f\x 
(l-a)^7«| s + a£A7$| s 



kM 



2nM 



kM 1 
12tt k r 



-5(f) , 



(l-a)d wl ( f\ s +aCA^>\ 



(5.2) 



Now, the static bulk modes are given by l|3.7|l : ^}(k l ,w) — e AI , y (fc l ) w 2 H^\ik\w\) and their w-derivatives on X 

are 01 : d w "fjfj\s = e^ u (k l ) ikC 2 (ikC). Equation (|5.2|) therefore gives the polarisations in terms of the brane 
stress-energy tensor as 



e 0Q (k l ) \{1 - a) h[ x) (ifc£) + iakCH 2 [1 \ikC) 
e tj [k l ) Ul-a) H[ X) (ikC) + i a kCH ( 2 1] (ikC) 
The bulk metric is then known and is : 



2 kM 1 

3 (27r)t ifc/: 2 
1 kM 1 
3(27r)iifc£ 2 



e 0l (fc l ) =0 
k j 



(5.3) 



7 M „(x\ w) = Tie 



d 3 k 



(27T) 



with 



7/*w (**)"' 



kM 



w 2 K 2 (k\w\) 



3£(27r)i kC [(1 - a)Kx(kC) + a kCK 2 (kC)} 



(5.4) 
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where cqq — 2, Co; = and Cj 



kikj/k 2 and where K u (z) is the modified Bessel function defined as K v {z) 



The hoQ Fourier component of the metric on the brane therefore reads 



hoo(k) = 7oo|s + 2— 
kM 



k 2 £{2n)i (1 + a) 



(1-fi) k£K (k£) 



3 (f + (5)^1 (fc£) + ak£K (k£) 



(5.5) 



and a similar expression for /iy . 

We then take the Fourier transform and integrate over angles (after elimination of the kikj term in hij by a suitable 
coordinate transformation). We then go to Schwarzschild coordinates by means of yet another infinitesimal coordinate 

transformation. Setting x = r / £ and recalling that we introduced the notation (^j-f^J ji = SttG, wc finally obtain 
the linearized metric on the brane, as created by a static point-like source, as : 



ds 2 | 4 = -(1 + 2U)dt 2 + (1 - 2V)dr 2 + r 2 (d0 2 + sin 2 



with 



and 



U 



GM 



l + —U & {x) 

07T 



V 



GM 



Ua{x) = (1 — a) / du sin(ux) 
Jo 

V a {x) = U & {x) -xd x Ua(x). 
Ua{x) is a well defined integral that we shall evaluate analytically. 



l + ^-V a (aO 



K {u) 



(5.6) 



(5.7) 



(1 + a)Ki(u) + auK Q (u) 



(5.8) 



VI. THE CORRECTIONS TO NEWTON'S LAW ON AN EINSTEIN BRANE FOR ALL DISTANCES 

When the bulk obeys the five dimensional Einstein equations, that is when a = 0, U a {x) in (|5.8|) reduces to 



U (x) = lim / du sm(ux)U (u)e~ eu with U (u) = ——7—7 

Ki(u) 



(6.1) 



By approximating Uq(u) by its asymptotic expressions, that is Uq(u) ~ 1 when u — * 00 and Uq(u) ~ — ulnu when 
u — > 0, it is a straightforward exercise [14J to see that lim^o Wo(x) = and that linix-^ % (x) = tt/2x 2 . Hence we 
recover the well-known result 0] (see also Q) that, at short distances the "correction" to the gravitational potentials 
diverges in £/r 



lim U = lim V = ~ — 

r/C^Q r/C^O 37T T Z 



(6.2) 



whereas at distances large compared with the characteristic scale £ of the anti-de Sitter bulk the correction is reduced 
by another £/r factor [lj, so that Q 



lim U = 



lim V = 



(6.3) 



However one can do better than just obtaining the asymptotic behaviours of the potentials. Indeed it turns out, 
see Fig. 1, that we can use for all u, to a good approximation 



Uo(u) « U (u) 



with 



Uq(u) = win 1 + - 



1 



Using that approximation, (|6.1|l can be integrated for all x to yield : 

x cos x — sin x cos x + x sin x 
Uq[x) w ci{x) -I si(x) 



7T 

2~? 



(6.4) 



(6.5) 



8 




Fig. 1. The behaviour of Uq(u) and the accuracy of the approximate analytic formula Uo(u) 
The solid curve is Uo(u) while the upper dotted line shows Uq/Uq. One sees that the 
relative error is small: 5; 2% at most. 



Log 10 [ ll Q (x) ] 
4 - 




Log 10 [x ; 



Fig. 2. A logarithmic plot of the function Uq(x). The short-dashed curve is tt/(2x 2 ) and the 

long-dashed curve is 1/x. 



and 



V (x) 



3a; cos x + (x 2 — 3) sin x 



x 

3x sin x 



(x 2 



ci(x) 

3) cos a; . , 
si{x) 



3tt 
2X 2 



1 

x 



(6.6) 



where si(a;) and ci(a;) are the integral sine and cosine functions, see [l4| . 

In Fig. 2, Uo(x) is plotted as a function of x — r/C. As one can see the transition from the 1/r behaviour at short 
distances to the 1/r 2 one at large distances is quite slow. This is due to the slow convergence of Uq{u) when u — * 0. 
As seen from the upper-most lines of Fig. 4 below, the Newtonian behaviour is recovered at r/C > 10. 
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VII. THE CORRECTIONS TO NEWTON'S LAW ON AN EINSTEIN GAUSS-BONNET BRANE 



When the bulk obeys the five dimensional Einstein Gauss-Bonnet equations, U a {x) in l|5.8[) is 

e + OO 



Ua(x) = (1 — a) j du sin(ux) U a (u) 

K (u) 



with Ua (u) 



auK (u) + (1 + a)K\{u) 
Again, by approximating Ua(u) by its asymptotic expressions, that is 

U a (u 

it is a simple exercise to see that 



(7.1) 



when u — ► oo , and U^Ou) ~ -u\nu when u — > 0, 

qu 1 + a 



7r 1 — a 



and lim U & {x) 



7r 1 — a 
2a; 2 1 + a 



Hence we recover the known result 0, that at distances large compared with the characteristic scale C of the anti-de 
Sitter bulk the gravitational potentials are corrected by a small 1 /r 2 factor, their precise expressions being 



GM 

lim U = 



lim V = 

r I C— >oo 



GM 



2 /l -a 



3 \1 + a 

1 — Q 



1 + a 



(7.2) 



When a = 0, (|7.2|l reduces to its "Einsteinian" expression ljti.311 . On the other hand the corrections remain finite when 
r/C — > so that the gravitational potentials keep their Newtonian behaviour : 



lim U = — - 
lim V = 



GM 


H 


-I 




r 






GM 


H 


l 


n-a.x 


r 




\ a )_ 



(7.3) 



contrarily to what happens in the Einstein case, see equation (|6.2() . 

Now, here again, one can do better by approximating U & {u) for all u. For example, for all < a < 1 (corresponding 
to the lower sign in l |2.15| l). a good approximation (to better than 1%) is 



U & (u) W U a (u) 



with 



U & (u) 



1 + aJ 1 + bu 



In 



(l + a) b 



(7.4) 



and b = 1 — y j+§ ■ Unfortunately we were not able to integrate (|7.1() analytically with that approximation. We hence 
used another approximate form for Ua (u) , valid for all a 7^ : 



Ua(u) ~ U a (u) with C/q 
where /? and 7 are constants such that 



l + a 



(7 2 ~l/3) 
(7-1/2) 2 



u In I 1 + - 

u 



l + a 



I3u 



u + /?7 



(7.5) 



= 



and (3 2 



l + a 



a (7-1/2) 



The relative error is found to be very small in both ranges of u 3> 1 and w«l for any a. At mildly small values of 
u, the error becomes as large as 25 %, but only for very large a (see Fig. 3 below). Using (|7.5II . I|7.1|l becomes 



U a {x) 



1 - a 
l + a 



(3x cos{(3x) — sin(/3x) .. /3a; sin(/3x) + cos(/3x) 

t ci(p x) + g si(p x) 



10 




Fig. 3. The behaviour of Ua(u) for several values of a and the accuracy of the approximate 
formula f/ a (li). The solid lines are (1 + a)Ua{u) for the cases a = 0.01, 0.3, 0.6, 0.99 and 
oo, from top to bottom. The dotted lines are the ratio Ua/Ua for the same values of a, 
also from top to bottom. The upper two dotted lines are almost degenerate in the figure. 



7T (3 

2x 2 x 



P 2 1 



ci(/3 7 x) sin(/3 7 x) — cos(/3 7 x) si(/3 7 x) 



(7.6) 



V a (a:) 



1 - a\ \ 2>(3x cos(x/3) + (x 2 /3 2 - 3) sm(/3x) , 

1 ' ^ '- ci(/3x) 



1 + a 



3(3x sin(Px) - {{? x 2 - 3) cos{(3 x) 3 it 3 j3 
si(/3x) + 7— 2 



+P 2 1 



i(/3 7 x) (shi(/3 7 x) — f3"f x cos(/3 7 a;) 
— si(/3 7 x) ( cos(/3 71) + (3 7 a; sin(/3 7 x 



2x 2 x 



(7.7) 



On the graphs below, see Fig. 4, we plotted U(r) and V(r) (as a function of x — r/£, and divided by the Newtonian 
potential Un = —GM/r). As one can see the gravitational potentials become closer and closer to their Newtonian 
values, for all distances, when a. — * 1. 



We therefore conclude that the Einstein "Lanczos Lovelock Gauss-Bonnet" equations in the bulk (rather than simply 
Einstein's) induce on the brane a better approximation to Newton's law and the usual linearized /owr-dimensional 
Einstein equations. This is all the more so as the coupling parameter a tends to 1. 

Let us be more precise : at short distances the potential U (r) remains in ljr when a ^ and hence C needs no 
longer be smaller than, say 200 ^m, the upper limit that present experiments 16] impose in the a = case for which 
U(r) cx 1/r 2 . 

As for the coupling constant G e = 2 j§ L G which appears in (7.3) we can identify it to Newton's constant, as measured 
in terrestrial laboratories (and known at present with a 0.15% accuracy, see e.g. |17]). 

If now, C is taken to be of geophysical size (say ~ 1 km —100 kms) then the transition from the short to large 
distance regimes does not contradict present experiments, whatever the value of a (not too close to zero), as it occurs 
"safely" in a region when gravity is poorly known. 

At large distances (r/C 3> 1) the potentials U(r) ~ V(r) are also in 1/r, with however a coupling constant G 
which differs from G e (identified with Newton's constant as measured in the lab). Now, there is no astvo-nomical 
evidence that G e should be equal to G, as only the product GM comes into play. There is on the other hand some 
astro-physical evidence that they should not be too different, otherwise solar models would go astray. If we tolerate 
that G and G e differ by, say, 10% then 0.85 < a < 1.15. 

To summarize : if £ is of geophysical size, and 0.85 < a < 1.15 (but =/= 1) then linearized gravity on an Einstein 
Gauss-Bonnet brane is compatible with Newtonian gravity and the light deflection experiments (which test the equality 
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U ~ V). 

We thus see that bulk Einstein Gauss-Bonnet gravity can relax the constraint on the value of C (that is, the 
characteristic size of the bulk) at a fairly resasonable price : a must be fine-tuned to be not too far from its critical 
value (1 - a) = (1 - 4a/C 2 ) = ^l + aA/3 = 0. 

Now, to see if such parameter ran ges can indeed be allowed, it will be necessary to analyze carefully their implications 
in the cosmological context as well |l8|. Investigations in this direction, however, are beyond the scope of the present 
paper. 
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